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An inelastic light (Raman) scattering study of nematicity and critical fluctuations in 
Ba(Fei_xCox) 2 As 2 (0 < x < 0.051) is presented. It is shown that the response from huctuations 
appears only in Big (x^ ~ y‘^) symmetry. The scattering amplitude increases towards the structural 
transition at Ts but vanishes only below the magnetic ordering transition at Tsdw < Tk , suggesting 
a magnetic origin of the huctuations. The theoretical analysis explains the selection rules and the 
temperature dependence of the huctuation response. These results make magnetism the favorite 
candidate for driving the series of transitions. 

PACS numbers: 74.70.Xa, 74.20.Mn, 74.25.nd, 74.40.-n 


I. INTRODUCTION 

Nematic huctuations and order play a prominent 
role in material classes such as the cuprates^, some 
ruthenates^ or the iron-based compounds^”^ and may be 
interrelated with superconductivity^”^^. In iron-based 
compounds^^’^^ signatures of nematicity have been ob¬ 
served in a variety of experiments, and the magneto- 
structural phase transition is among the most thoroughly 
studied phenomena. When Fe is substituted by Go in 
BaFe 2 As 2 the structural transformation at precedes 
the magnetic ordering at Tsdw < Tg [Ref. 14]. The 
nematic phase between Tg and Tsdw is characterized 
by broken C 4 symmetry but preserved 0(3) spin ro¬ 
tational symmetry (no magnetic order). Nematic huc¬ 
tuations are present even above Tg in the tetragonal 
phase as has been demonstrated in studies of the elas¬ 
tic constants^^. In strained samples, one observes orbital 
ordering in the photoemission spectra^^ and electronic 
nematicity by transport^’However, the fundamental 
question as to the relevance of the related spin^^, charge^^ 
or orbitaF’^^’^^ huctuations remains open. In fact, it is 
rather difficult to derive the dynamics and momentum 
dependence of the critical huctuations with hnite char¬ 
acteristic wavelengths^and to identify which of the 
ordering phenomena drives the instabilities. 

Raman scattering provides experimental access to all 
types of dynamic nematicity but only the charge sector 
has been studied in more detaiF^’^^’^^. However, also in 
the case of spin-driven nematic order the technique can 
play a prominent role for coupling to a two-spin opera¬ 
tor whereas a four-spin correlation function is the lowest 
order contribution to the neutron cross section^. We ex¬ 
ploit this advantage here and study the low-energy Ra¬ 
man response of Ba(Fei_xCox) 2 As 2 experimentally and 


interpret the results in terms of a microscopic model for 
a spin-driven nematic phase. In addition to the tempera¬ 
ture dependence^^ we address the spectral shape and the 
selection rules enabling us to explain the structural and 
magnetic transitions in a unified microscopic picture. 

We study Ba(Fei_xCox) 2 As 2 single crystals having 
X = 0, X = 0.025, and x = 0.051 as a function of photon 
polarization in the temperature range 4.2 < T < 300 K. 
For the symmetry assignment we use the 1 Fe unit cell 
making the fluctuations to appear in Big symmetry. We 
use the appearance of twin boundaries and of the As Aig 
{x^ y‘^) phonon line as internal thermometers for the 

structural and the magnetic phase transitions, respec¬ 
tively. In this way, Tg and Tsdw can be determined with 
a precision of typically ±0.2 and ±1K, respectively. 


II. EXPERIMENT 

The single crystals of undoped and Go-substituted 
Ba(Fei_xCox) 2 As 2 were grown using a self-flux technique 
and have been characterized elsewhere^^. The cobalt con¬ 
centration was determined by microprobe analysis. Tg 
and Tsdw a.re close to 134 K in the undoped sample and 
cannot be distinguished. At nominally x = 0.025 we 
find Tg = 102.8 ± 0.1 K and Tsdw = 98 ± IK by di¬ 
rectly observing the appearance of twin boundaries and 
a symmetry-forbidden phonon line, respectively (see Ap¬ 
pendix A for details). The extremely sharp transition at 
Tg having ATg ^ 0.2K indicates that the sample is very 
homogeneous in the area of the laser spot. 

The experiments were performed with standard light 
scattering equipment. For excitation either a solid state 
laser (Goherent, Sapphire SF 532-155 GW) or an Ar ion 
laser (Goherent, Innova 300) was used emitting at 532 
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FIG. 1. Raman response Rx"(Q,T) (raw data after division by the Bose-Einstein factor) of Ba(Feo. 975 Coo.o 25 ) 2 As 2 in (a) Big 
above and (b) below Tg and (c) Aig symmetry at temperatures as indicated. The initial slopes shown in (a) and (c) as grey 
arrows are proportional to the static two-particle lifetime in symmetry // = Aig, Big. (d) Raman relaxation rates ro,^(T) in 
Aig (blue circles) and Big (red diamonds) symmetry as a function of temperature. The fluctuation range Tg < T < Tf and 
the nematic phase Tsdw < T < Tg are indicated in green and magenta, respectively. 


or 514.5 nm, respectively. The samples were mounted on 
the cold finger of a He-flow cryostat in a cryogenically 
pumped vacuum. The laser-induced heating was deter¬ 
mined experimentally (see Appendix A) to be close to 
1 K per mW absorbed power. The spectra represent the 
response Rx'^{^^T) (/i = Big^ A 2 g and B 2 g) that is 
obtained by dividing the measured (symmetry resolved) 
spectra by the Bose-Einstein factor {1 + n(T, 11)} = 
[1 —exp(—x|^(ll,T) is the imaginary part of 
the response function, and R is an experimental constant 
that connects the observed photon count rates with the 
cross-section and the van Hove function and accounts for 
units. For simplicity the symmetry index fi is dropped in 
most of the cases. The symmetry selection rules refer to 
the 1 Fe unit cell (see insert of Fig. 6 (b) in Appendix A) 
which is more appropriate for electronic and spin excita¬ 
tions. 


III. EXPERIMENTAL RESULTS 

Fig. 1 shows the Raman response Rx"(ll,T) for 
Ba(Feo. 975 Coo.o 25 ) 2 As 2 for various temperatures in Aig 
and Big (IFe per unit cell) symmetry. B 2 g spectra were 
measured only at a few temperatures and found to be 
nearly temperature independent in agreement with pre¬ 
vious data^^. Results for other doping levels x are shown 
in Appendix B. The spectra comprise a superposition 
of several types of excitations including narrow phonon 
lines and slowly varying continua arising from electron- 


hole (e-h) pairs; hence the continuum reflects the dy¬ 
namical two-particle behavior. The Aig and Big spectra 
predominantly weigh out contributions from the central 
hole bands and the electron bands, respectively^^’^^. The 
symmetry-dependent initial slope ro,^(T) (/i = Aig, Big, 
B 2 g) [see Fig. 1 (a) and (c)] can be compared to trans¬ 
port data. [ro,^(T)]“^ corresponds to the static transport 
relaxation rate Fo,^(T) of the conduction electrons^^”^^. 
The memory function method facilitates the quantitative 
determination of the dynamic relaxation F(U,T) in ab¬ 
solute energy units^^. The static limit can be obtained 
by extrapolation, Fo,^(T) = F^(U —0,T) [see Appendix 
C]. In Fig. 1 (d) we show the result for x = 0.025 corre¬ 
sponding to the spectra of Fig. 1 (a), (b) and (c). The 
results for all doping levels studied are compiled in Fig. 9 
in Appendix C and compared to the scattering rates de¬ 
rived from the resistivities^^. 

Fig. 1 (d) displays one of the central results: Above ap¬ 
proximately 200 K Fo,^(T) varies slowly and similarly in 
both symmetries. The more rapid decrease of To^Big{T) 
below 200 K is accompanied by a strong intensity gain in 
the range 20-200 cm“^ [see Fig. 1 (a)] as observed before 
in similar samples^^’^^. The intensity gain indicates that 
there is an additional contribution superposed on the e- 
h continuum which, as will be shown below, arises from 
fluctuations. Therefore, the kink in To^BigiT) is labeled 
Tf and marks the crossover temperature below which ne¬ 
matic fluctuations can be observed by Raman scattering. 
At least for low doping, Tf is relatively well defined. The 
kink allows us to separate the two regimes of the low- 
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FIG. 2. Fluctuation contribution to the Raman response of 
Ba(Feo. 975 Coo.o 25 ) 2 As 2 (a) above Tg and (b) between Tg and 
Fsdw- The red lines are theoretical predictions on the basis 
of Aslamazov-Larkin diagrams^^ describing the exchange of a 
pair of fluctuations (for details see Appendix D). 


energy response above and below Tj as being dominated 
by carrier excitations and fluctuations, respectively. 

The additional Big signal below Tf has to be treated 
in a way different from that in Aig symmetry and in 
Big above Tf. Since it is rather strong it can be sepa¬ 
rated out with little uncertainty by subtracting the e-h 
continuum. We approximate the continuum at Tj by an 
analytic function which is then determined for each tem¬ 
perature according to the variation of the resistivity and 
the Aig spectra and subtracted from all spectra at lower 
temperatures. The details are explained in Appendix E. 
The results of the subtraction procedure are shown in 
Fig. 2. The response increases rapidly towards Tg with¬ 
out however diverging, and the maximum moves to lower 
energies. 


IV. THEORY 


We first compare the data to the theoretical model for 
thermally driven spin fluctuations associated with the 
striped magnetic phase ordering along = (tt, 0) or 
= (0,7r). In leading order two noninteracting fluc¬ 
tuations carrying momenta Q and — Q are exchanged. 
Electronic loops (see Eigs. 3 and 10) connect the photons 
and the fluctuations and entail Q-dependent selection 
rules that were derived along with the response 
in Ref. 22 and are summarized in Appendix D. In brief, 
since the response results from a sum over all electronic 
momenta close to the Eermi surface cancellation effects 
may occur if Q connects parts on different Eermi surface 
sheets having form factors 7 ^(k) with opposite sign. Eor 
the ordering vectors (tt, 0) and (0, tt) the resulting selec¬ 
tion rules explain the enhancement of the signal in Big 
symmetry and its absence in the Aig and B 2 g channels. 
In contrast, for ferro-orbital ordering with Q = (0, 0) as 
found in EeSe^ the fluctuation response would appear in 
all symmetries. 

However, the lowest-order diagrams alone can only ac¬ 
count for the spectral shape whereas the variation of the 
intensity around Tg remains unexplained. In order to 
describe this aspect, we consider the interaction of fluc¬ 
tuations among themselves and with the lattice, all of 
which becomes crucial in the treatment of spin-driven 
nematicity^^. 

The interactions between spin fluctuations can be rep¬ 
resented by a series of quaternion paramagnetic couplings 
mediated by fermions inserted into the leading order 
Aslamazov-Larkin diagrams as shown in Eig. 3. The in¬ 
serted fermionic boxes effectively resemble the dynamic 
nematic coupling constant g of the theory. 

We have analyzed the problem by extending SU (2) ^ 
SU{N) and taking the large N limit. Eor small frequen¬ 
cies Cl and in the large-V limit, after re-summing an infi¬ 
nite number of such box-like Aslamazov-Larkin diagrams, 
the Raman response function Rsigi^) reads, 

RBlAm = Ro,BlAm [1 + 5X^'em(0)] • (1) 

Eq. (1) states that the Raman response is proportional 
to the electronic contribution to the susceptibility of the 
nematic order parameter. 
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As a surprise, the fluctuations do not disappear di¬ 
rectly below Tg [Eig. 2 (b)] as one would expect if long- 
ranged order would be established. Rather, the intensity 
decreases continuously and the maximum stays approxi¬ 
mately pinned implying that the correlation length does 
not change substantially between the two transitions at 
Tg = 102.8 ±0.2 K and Tsdw = 98T1K. The persistence 
of the fluctuations down to Tsdw strongly favors their 
magnetic origin. 


Xmag(^) represents the magnetic susceptibility that di¬ 
verges at Tsdw- For g ^ 0 Xnlm(O) Fas a Curie-like 
|T _ T*|-i divergence at T* > Tsdw- 
If the spins (or charges) couple to the lattice the 
susceptibility of the nematic order parameter is given 
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where Agi denotes the magneto-elastic coupling, and Cq is 
the bare elastic constant. Obviously, Xnem(O) diverges at 
higher temperature than Xnem(O)* We identify > T* 
with the structural transition and conclude that the Ra¬ 
man response (Eq. 1) develops only a maximum rather 
than a divergence at Tg in agreement with the experi¬ 
ment. 

Close to T 5 , we expect Eq. ( 1 ) to hold qualitatively also 
inside the nematic phase, Tsdw < T < Tg. We argue^^ 
that Xnem(O) accordiug to Eq. (1), the Raman ampli¬ 
tude is smaller than in the disordered (tetragonal) state 
but different from zero. This explains the continuous re¬ 
duction of the Raman response of spin fluctuations upon 
entering the nematic state. One can also show that the 
Aig response gets even further suppressed if one includes 
collisions between the fluctuations. 



FIG. 3. Higher order Aslamazov-Larkin diagrams for interact¬ 
ing fluctuations. The paramagnetic coupling between fluctu¬ 
ations mediated by fermions is obtained by inserting quater¬ 
nions into the leading order diagram. The re-summed box 
Ba /3 is shaded grey. The first index of the matrix B denotes 
the type a — XjY of entering spin fluctuations, and the sec¬ 
ond index the type of exiting spin fluctuations. 


V. DISCUSSION 

As shown in Eq. ( 1 ) the full Raman response is 

proportional to the bare response and to the elec¬ 

tronic nematic susceptibility Xnlm(O)- Hence, the spec¬ 
tral shape is essentially given by that is there¬ 

fore used in Eig. 2 to fit the data, whereas the intensity is 
dominated by the prefactor |T—T* Since the theoret¬ 
ical model is valid only in the limit of small frequencies we 
argue that the initial slope reflects the temperature de¬ 
pendence of the intensity and is proportional to XnemCO)^ 
at least close to the transition. Eor generally reflecting 
the spectral shape above Tsdw (Eq. 1), Ro^Bigi^^T) en¬ 
ables us to directly extract the initial slope of the exper¬ 
imental spectra by plotting Ro,bi ^(^5 T)/U for all tem¬ 
peratures (see Appendix E). These results are compiled 
in Eig. 4 along with the variation of Xnlm (05 expected 
from mean-field theory. Eor low doping, we find qualita¬ 
tive agreement in the ranges Tsdw < T < Tg and Tg < T. 
Eor higher doping the interactions between fluctuations 
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FIG. 4. Temperature dependence of the initial slope of the 
fluctuation response. The initial slope is given in arbitrary 
units since only the temperature dependence matters. The 
error bars represent the differences originating in the sub¬ 
traction of the e-h continuum. The grey curves represent 
|T _ T*|-i with T* < Ts (see text). 


become dominant and the mean field prediction breaks 
down [Fig. 4(c)]. 

The fluctuations were also studied at various other 
doping levels in the range 0 < x < 0.085. Up to 6 . 1 % Co 
substitution fluctuations were observed. In contrast to 
other publications^^ we were not able to clearly identify 
and isolate the response from fluctuations at 8.5%. The 
results up to 5.1% are unambiguous and are compiled in 
Fig. 5. The fluctuations can be observed over a tempera¬ 
ture range of approximately 70-100 K. This is more than 
in most of the other experiments on unstrained samples 
and comparable to what is found in the cuprates^^’^^. 


VI. CONCLUSIONS 

The detailed experimental and theoretical study of the 
light scattering response in Ba(Fei_xCox) 2 As 2 reveals a 
broad range of spin fluctuations obeying Big selection 
rules. The selection rules can be explained only for a fi¬ 
nite ordering vector Q = ( 7 r, 0 ). Any type of order with 
Q = (tt, tt) or Q = (0, 0) such as ferro-orbital order is not 
compatible with the experiment. By observing the twin 
boundaries and the As phonon intensity we are able to 
determine the structural and magnetic transition temper¬ 
atures with unprecedented precision. This observation 
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FIG. 5. Phase diagram of Ba(Fei_xCox) 2 As 2 . The full lines 
limiting the nematic phase (magenta) and the blue squares 
representing the transition temperature Tc of superconducting 
samples were derived in Ref. 14. Grey diamonds represent 
doping and temperature positions of the current Raman data. 
The red-green field between Ts and Tf represents the initial 
slope of the spectra according to the color scale on the right. 


allows us to conclude that the intensity of the fluctua¬ 
tion response is maximal at Tg and vanishes at Tsdw- 
The divergence of the intensity expected at Tg from the 
electronic nematic susceptibility alone is shifted to lower 
temperature due to magneto-elastic coupling. Therefore, 
only an intensity maximum is observed at Tg. This fact 
together with the observation that the signal disappears 
at Tsdw supports the spin-driven nematic phase scenario. 
Hence magnetism is likely to be behind the transitions at 
least in Ba(Fei_xCox) 2 As 2 and makes its fluctuations a 
candidate for driving superconductivity. 
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FIG. 6. Determination of the spot temperature, (al) and 
(a2) show the sample surface and a zoom in thereof. In (a2) 
the approximate laser spot is indicated schematically. (a3) 
shows the crystallographic axes of the tetragonal phase, (b) 
Ag phonon around Tsdw • Above Tsdw the symmetry leakage 
of the phonon is negligibly small. Only below Tsdw the inten¬ 
sity becomes appreciable, (c) Image of the illuminated spot as 
a function of temperature. The horizontal lines in (cl)-(c3) 
result from twin boundaries. In addition to the appearance 
of twin boundaries the reversible adsorption of residual gas 
atoms and molecules starts instantaneously and enhances the 
stray-light in the spot. 


Appendix A: Determination of the spot temperature 

In Figs. 2, 7 and 8 we show that the response from 
fluctuations is maximal at Tg and then decreases. For 
X = 0 the decrease is very rapid, at x = 0.025 and 
0.051 the fluctuations disappear only below Tsdw- Since 
AT = Tg— Tsdw is small close to zero doping, the laser- 
induced heating has to be determined precisely. In ad¬ 
dition, a large temperature gradient in the spot would 
lead to a substantial reduction of the maximal fluctua¬ 
tion intensity. Great care was therefore taken to keep 
the temperature gradient in the spot small and to de¬ 
termine the spot temperature and to calibrate it against 
intrinsic thermometers. The calibration is possible since 
twin boundaries develop below Tg in the samples with 
X = 0.025 and 0.051 facilitating a very precise determi¬ 
nation of Tg. First we studied the effect of increasing 
laser power Pl at different holder temperatures Th on 
the twin pattern that can be seen, e.g., in Fig. 6 (cl). 
In this way the laser heating ATl was determined to be 
1 ± O.lK/mW for a spot diameter d = 50/im. (Note 
that ATl scales as d~^ and not as d~‘^.) Next we heat 
the sample slowly through Tg using Pl = 0.3 mW as 
shown in a series of snapshots in Fig. 6(cl)-(c6). The 
twin boundaries appear as horizontal lines and are most 
pronounced in (cl). With increasing temperature they 
“melt” and finally disappear completely at 102.9 K (ex¬ 
trapolated sample temperature for Pl = 0), and we iden¬ 
tify Tg = 102.9 K. 

For estimating Tsdw we analyze the phonons. The Aig 
As vibration was reported to appear in B 2 g symmetry be¬ 
low Tg^^ . (We maintain the tetragonal 2 Fe unit cell here 
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FIG. 7. Raman response Rx"{^,T) (raw data) of BaFe 2 As 2 
in (a), (b) Big and (c) Aig symmetry above and below the 
structural transition Ts at temperatures as indicated. 


FIG. 8. Raman response Rx"{^,T) (raw data) of 
Ba(Feo.949Goo.o5i)2As2 in (a), (b) Big and (c) Aig symmetry 
above and below the structural transition Ts at temperatures 
as indicated. The inset in (b) shows that the SDW gap starts 
opening within 5 K below Tsdw • 


as opposed to the main text to avoid confusion with the 
usual phonon assignment. In the proper orthorhombic 
4Fe unit cell applying below Tsdw the phonon switches 
to Ag symmetry, and B 2 g symmetry is not accessible any 
further with in-plan polarizations.) Our precise temper¬ 
ature determination shows for x = 0.025 that the anoma¬ 
lous intensity does not appear at Tg. Rather the phonon 
anomaly appears only at approximately 97 K as shown in 
Fig. 6 (b). According to the phase diagram the magnetic 
transition is offset by approximately 4-5 K at x = 0.025. 
This is actually not unexpected for a phonon that is 
not coupled to the lattice distortion by symmetry^^. By 
measuring the B 2 g intensity of the Aig phonon we can 
therefore identify the magnetic transition temperature 
and find Tsdw = 98 ± 1 K. 

For X = 0.051 we find Tg = 61.0 ± 0.2 K and Tsdw = 
51 ± 2 K. Here, the Aig phonon appears already above 
Tsdw 5 and we identify Tsdw with the strongest increase 
of the intensity. In addition, we know the width of the 
nematic phase from the phase diagram^^ (Fig. 5) and 
find an anomaly of Fo,^(T) close to Tsdw [see Fig. 9 (d)]. 
Hence the relevant temperatures are known with high 
precision. 


Appendix B: Results at x = 0 and x — 0.051 


Figs. 7 and 8 show the experimental results for x = 0 
and x = 0.051. Atx = 0 the two transitions Tsdw 
and Tg either coincide or are too close to be observed 
separately while the response of the SDW phase can be 
identified clearly as observed earlier^^’^^. At x = 0.051 
the fluctuations can be separated out in the usual way 
as described below. If an extraction is attempted in a 
similar way at x = 0.085 the variation with tempera¬ 
ture cannot be described with Aslamazov-Larkin-type of 
fluctuations. Although the response increases slightly to¬ 
wards lower temperature^^ and the elastic constants may 
still indicate an instability up to 9% Co substitution^^ we 
do not feel comfortable to extract parameters in this case. 
The results for Fo,^(T) are compiled in Fig. 9. 
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FIG. 9. Static Raman relaxation rates ro(T) of (a) BaFe 2 As 2 , 
(b) Ba(Feo.975Coo.o25)2As2, and (c) Ba(Feo.949Coo.o5i)2As2. 
Fo(T) is derived formally via the memory function method^® 
as described in appendix C. Above the onset temperature 
of the fluctuations Tf the results in both Aig (blue) and 
Big (red) symmetry return results similar to those from the 
resistivity^^ (right ordinate). Using a Drude model, the resis¬ 
tivities p(T) can be converted into scattering rates. At T/ the 
temperature dependence in the Big symmetry becomes much 
stronger, (d) If the resolution in temperatures is very high, 
one finds anomalies of Fo,^(T) at Ts and Tsdw which facili¬ 
tates the independent determination of Ts and Tsdw directly 
from the electronic Raman spectra. 



FIG. 10. Scattering from fluctuations, (a) Examples of 
Aslamazov-Larkin diagrams describing light scattering from 
critical fluctuations with momentum zbqc and energy ujm- 
The dashed lines represent the magnetic fluctuations, the 
full lines the fermionic propagators, (b) Hole- (grey) and 
back-folded electron-like (black) Fermi surfaces intersecting in 
the hot-spots (yellow circles). The selection rules can be de¬ 
duced by considering cancelation effects arising from different 
hot-spot contributions inside the fermionic loops described in 
Eq. (Dl). The first and second row of (c), (d) and (e) show 
the signs and nodes of the first and second order Aig, Big, and 
B 2 g Brillouin zone harmonics that indicate where cancellation 
effects can and cannot be expected. The qc vectors for (tt, 0) 
and equivalent fluctuations are indicated by full and broken 
arrows, respectively. The last row shows the vertices derived 
from the second derivative of tight-binding band structure (ef- 
fectice mass approximation) of Graser et al. (Ref. 37). These 
vertices provide the best estimate for the sensitivity on the 
Eermi surface^^. The Aig vertices for the hole and the elec¬ 
tron bands are predominantly negative (blue) and positive 
(red), respectively. The effective mass approximation shows 
that the Aig response will be dominated by the second or¬ 
der vertex cos kx cos ky rather than the lowest order one as 
already pointed out in Ref. 28. 








































Appendix C: Memory function and static relaxation 

rates 


In Fig. 1 (d) symmetry-dependent static relaxation 
rates ro,^(T) are shown for /i = Aig and Big^ 


h 

roAT) 


ro,M(T) 


dRx''{n,T) 

dn 


-1 


(Cl) 


0=0 


Since the overall intensity of the spectra is not known 
in absolute units the experimental constant i?, to which 
the initial slope ro,^(T) is proportional, cannot be pinned 
down. Therefore, one needs additional information if 
one is interested in energy units for ro,^(T). Only then 
ro,^(F) can be compared to transport data. This prob¬ 
lem was solved a while ago by adopting the memory func¬ 
tion method^^’^^ for Raman scattering^^ . Then ro,^(T) 
can be derived by extrapolating the dynamic Raman re¬ 
laxation rates r^(^l,T) = The results for all 

doping levels are compiled in Fig. 9. 

If a Drude model is applied the resistivities p{T) can 
be converted into static scattering rates. Using a plasma 
frequency close to 1 eV in rough agreement with optical 
data"^^, the analysis shows that the Raman and trans¬ 
port results are compatible above a doping dependent 
temperature Tf that is identified here with the onset of 
fluctuations in agreement with results from other meth¬ 
ods. Transport and Raman scattering agree to within 
the experimental precision, possibly indicating the com¬ 
mon origin of the electronic relaxation on the electron 
and hole bands. 


Appendix D: Aslamazov-Larkin Diagrams and 
Selection Rules 

The coupling of visible light to critical fluctuations 
with wavevectors |qc| = Q > 0 and energy (mass) uj^n is 
possible only via the creation of two excitations with op¬ 
posite momenta warranting zero net momentum transfer 
applying for photon energies in the eV range [Fig. 10 (a)]. 
This process can be described by Aslamazov-Larkin (AL) 
diagrams^^. We assume a simplified model of the Fermi 
surface. The central sheet is a circular hole-like pocket 
around the F point [grey circle in Fig. 10(b)]. The two 
electron-like elliptical pockets with the principle axes ro¬ 
tated by 90^ are centered at the X (Ttt, 0) and Y (0, Ttt) 
points of the 1 Fe BZ. If they are backfolded they inter¬ 
sect with the central hole band as indicated by yellow 
circle in [Fig. 10(b)]. The fluctuation contribution to 
the Raman spectrum has been analyzed by Caprara and 
coworkers for the cuprates^^ and arises from the AL di¬ 
agrams shown in Fig 10(a). The selection rules can be 
deduced by considering cancelation effects arising from 
different hot-spots within the fermionic loop as shown 
in Fig. 10 (a). Even if the entire Fermi surface is taken 
into account the selection rules still work in the Fe-based 


materials. For instance, in either case full cancellation is 
found for B 2 g symmetry^^. 

Explicitly written out, the fermionic loop is given 

by22-24,41 


^i,/x(qc5 ^m) — i^ci ^9 ^m) T ^m)i 

e^%,,n,com)=Tj2 [ l^Gr{k,Sn-n)Gr{k,en) 

n 

=T^ /* - n)Gi(k,£„) 

n 

X G^r (h qc9 U T (^^) 


( 2 ). 


where 7 ^ is the form factor (/i = Big^ Aig etc.), and 
Gi is the electron propagator on band i = F, X, Y. 
is the electronic energy and Q is the energy difference 
between the incoming and scattered photons. Experi¬ 
mentally, pure symmetries can be obtained from linear 
combinations of the response measured at appropriate 
polarizations of the incoming and scattered photons 
and 63 . 

Eor illustration purposes the fermionic loop 0 is ap¬ 
proximated in the hot-spot approximation. Hot-spots 
are regions in momentum space where both k and k ± qc 
lie on the Eermi surface [Eig 10(b)]. Since the loop 0 
contains the symmetry factor 7 (k) linearly inside the 
momentum integral the sign of 7 (k) is crucial. If 7 (k) 
changes sign for different hot spots connected by qc 
(Eig. 10(c), (d), and (e) for Ai^, Big^ ^ 2 ^, respectively) 
there will be full or partial cancelation within 0. Eull 
cancelation is observed for the first two (and also higher) 
orders of B 2 g symmetry [Eig. 10(e)]. In contrast, 7 (k) 
does not change sign across different hot-spots for the Big 
channel. Consequently, in Big and B 2 g the fluctuations 
are Raman active and inactive, respectively. 

The Aig symmetry is more complicated in that the 
first order contribution, proportional to co^(kx) Y co^{ky) 
[upper row of Eig. 10(c)], would be as strong as the 
Big contribution [Eig. 10 (d)] whereas the second order 
contribution (co^{kx) co^{ky)) [second row of Eig. 10(c)] 
shows cancelation. Eor clarifying the relative magnitude 
of the two orders we analyze the effective mass vertices 
on the Eermi surfaces (second derivative or curvature of 
the band structure), that are the best approximations for 
the sensitivity away from resonances, in a way similar to 
what was proposed in Ref. 28. The last row of Eig. 10 (c) 
shows that the band curvatures corresponding to the Aig 
vertex 


7i,/llg(k) 


dkxdkx dkydky 


(D2) 


on the Eermi surface of the hole and the electron bands (i) 
are predomininantly negative and positive, respectively, 
as expected already for simple parabolic bands with 
masses rrih ~ —rrie although there are various near nodes 
on both bands. This result shows that cos{kx) cos{ky) is 
the leading order. We note that cos{kx) cos{ky) predicts 







9 


c 

sS 

p 

d 

q: 


0.6 

0.4 

0.2 

0.0 


(a) ' 

' 1 1 _ 

/ 

x = 0.025 _ 

/ 

- 300 K 

- If 

200 K - 

/ 'i r 

150 K 

t _1_ 

- e-h fit 

_I_I_ 


0 200 400 600 800 

Raman Shift Q (cm'^) 


- 125 

K 

- 115 

K 

- 105 

K 

- 103 

K 

- 102 

K 

- 101 

K 

CD 

CD 

K 



Raman Shift Q (cm'^) 




FIG. 11. Determination of the e-h continuum and subtraction from the full response in Ba(Fei_xCox) 2 As 2 {x = 0.025). (a) 
The analytical function (red dashes) provides a reasonable ht at 300 K. At lower temperatures the fluctuations emerge above 
the continuum and the analytical functions lie below the data (yellow and green dashes), (b) The analytical model is varied 
so as to reproduce the temperature dependence of the resistivity as shown as orange diamonds in panel (d). (c) Dynamical 
relaxation rates F(Q,T) derived from the synthetic spectra in panel (b). The zero-energy extrapolation values of Fc(0,T) are 
plotted as orange diamonds in (d). The Aig and Big data are taken from Fig. 1. 


a stronger mixing of the particle-hole response from the 
electron and hole bands than cos{kx)cos{ky) as already 
outlined by Mazin et al Ref. 28. 


Appendix E: Subtraction of the continuum 


The fluctuation response is superposed on the particle- 
hole continuum that essentially reflects symmetry- 
resolved transport properties'^. Since the contribution 
of the fluctuations is relatively strong here they can be 
isolated with little uncertainty. The simplest way is to 
use the continuum at or slightly above the crossover 
temperature Tf and subtract it from all spectra mea¬ 
sured below Tf. This was sufficient for ErTes^^ but cre¬ 
ated negative intensities in the case of La 2 -xSrxCu 04 ^^. 
Here, we wish to compare the temperature dependence 
of the fluctuations to a theoretical prediction and have 
to improve on the subtraction of the continuum. To this 
end we make the analytical phenomenology for the Big 
continuum temperature dependent in a way that yields 
(X p{T). This seems sensible since the propor¬ 
tionality holds for the Aig results in the entire temper¬ 
ature range above Tsdw and for the Big spectra above 
Tf. Fig. 11 shows the steps and checks necessary for the 


procedure. The analytical function used reads 

Xcont(^5^) = [<^i + <^2 • T] tanh + 

which obeys x"o„t{-n,T) = -Xconti^^B as required by 
causality, o^i, 0 ^ 2 , Pi and P 2 depend only on doping x. 
For X = 0.025 we used ai = 0.82379, 0^2 = —0.00138, 
Pi = —0.00923, and P 2 = 0.00028. fo(^) is a fitting pa¬ 
rameter that is selected in a way that the inverse slope 
Fc(0,T) of x^q^^( 11,T) follows the resistivity (orange di¬ 
amonds in Fig. lid). If a constant continuum is used 
the fluctuations can be isolated in a qualitatively similar 
fashion. However, the experimental data in Fig. 2 vary 
more slowly close to T^. 

Below Ts the uncertainties increase since surface layers 
accumulate rapidly in the presence of twin boundaries 
where the surface assumes a more polar character. This 
can be seen directly in Fig. 6 (c). 


Appendix F: Initial slope 

For being a causal function the Raman response is an¬ 
tisymmetric and, as long as there is no gap, linear around 
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FIG. 12. Initial slope of the Raman response of 
Ba(Fei-xCox) 2 As 2 {x = 0.025) below T/. This hgure is a 
reproduction of Fig. 2 with the fluctuation response divided 
by the energy Q, plotted against a logarithmic energy scale. 


the origin. Then Eq. (Cl) can be approximated as 


roAT) 


dRxm,T) 


lim 

o—^0 


on 

n 


0=0 


(FI) 


The temperature dependence (not the magnitude) of the 
initial slope can then directly be read off a graph if the 
response is divided by the energy ft and plotted against 
a logarithmic energy scale. 

If R was known ro(T) could be determined directly. 
With R unknown only the relative change can be derived 
in this way. Fig. 12 shows that the fits reproduce the 
overall data rather well at low energy. The phenomeno¬ 
logical curves can be extended to arbitrarily low energies 
providing a simple way to directly visualize the temper¬ 
ature dependence of to(T). Fig. 12 shows also that the 
experimental data close to zero energy are not very sta¬ 
ble. This problem arises from accumulating surface layers 
and the influence of the laser line. Therefore, the error 
bars become excessively large if the slope is directly ex¬ 
tracted from the data. Here we use a wide spectral range 
to improve the reproducibility. 
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